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ABSTRACT 

Ran ©arch  performed  by  Robert  Davison 
Under  The  Supervision  of  Dr#  rt#  J#  McHichols 

Thin  paper  develops  a  method  to  allocate  reliability 
to  system  elements  where  a  system  reliability  requirement 
must  be  met  or  a  total  coat  constraint  is  imposed#  Parallel 
and  serially  arranged  elements  are  considered  with  each 
element's  probability  density  function  of  time  to  failure 
being  approximated  by  the  two  parameter  Woibull  distribution 
The  problem  is  approached  as  a  minimization  of  total  system 
cost  or  system  unreliability  and  Lagrange  multipliers  prove 
to  be  useful  as  a  solution  technique#  A  cost  function  is 
developed  to  relate  the  cost  of  and  element  to  the  relia¬ 
bility  level  achieved  in  that  element.  Both  the  element's 
cost  function  and  reliability  expression  must  be  continuous 
variables#  numerical  examples  are  shown  for  each  allocation 
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CHAPT33  I 


INTRODUCTION  TO  RELIABILITY  ALLOC  ATT  ON 

During  several  phases  of  design  and  development  of  sys- 
ten3  for  Army  use,  it  is  necessary  to  translate  the  system’s 
and  subsystems’  availability,  maintainability  and  readiness 
requirements  into  quantitative  reliability  requirements  for 
the  many  system  sublevels.  It  is  mandatory  that  both  real¬ 
istic  and  consistent  reliability  specifications  be  assigned 
to  all  the  components  and  or  subsystems  that  comprise  the 
desired  operational  package.  The  process  of  this  assignment 
is  most  commonly  termed  reliability  apportionment  or  relia¬ 
bility  allocation. 

The  prime  factor  to  be  considered  in  any  reliability 
allocation  is  that  of  assuring  the  total  system  reliability 
requirement  is  met  once  the  component  reliabilitios  are  ob¬ 
tained  and  the  system  assembled  for  operation.  More  fre¬ 
quently  than  not  it  is  a  cumbersome  if  not  impossible  task 
to  exactly  express  the  reliability  of  the  system  as  a  func¬ 
tion  of  the  components'  and  subsystems'  reliabilities  in  a 
nathmatical  statement.  Several  factors  contribute  to  this 
difficulty.  Early  in  the  design  phase  3one  subsystem  de¬ 
signs  are  not  nfroaon;t  with  respect  to  configuration  and 
complexity.  Essentiality  of  components  to  system  success  is 
not  well  known.  Certain  subsystems  are  O’Jy  needod  intor- 
m.ifcfcantly  or  for  a  portion  of  the  mission  while  others 
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might  bo  required  for  the  duration.  The  uoago  environment 
is  not  always  well  defined  and  indeed  the  syaton  night  im¬ 
pose  additional  stress  on  its  components  os  a  consoquonoe  of 
its  operation.  Where  a  system  utilises  state-of-the-art 
components,  the  foregoing  problems  are  compounded  by  the 
lack  of  historical  data  concerning  reliability.  To  make  the 
reliability  allocation  even  more  complicated  is  one  addition¬ 
al  constraint;  economy.  Today,  with  the  financial  scrutiny 
being  exercised  in  all  categories  of  government  expenditure, 
it  is  important  that  the  most  effective  system  possible  be 
obtained  for  the  resources  available. 

In  the  past  two  decades  much  work  has  gone  into  the 
development  of  methods  to  allocate  reliability  to  system 
sublevels.  Earlier  allocation  methods  concerned  themselves 
solely  with  the  problem  of  satisfying  a  reliability  or  fail¬ 
ure  rate  specification.  ‘To  specific  consideration  for  coat 
was  reflected  in  the  models  used  to  allooate  system  relia¬ 
bility  or  failure  rate.  Although  design  specifications 
:  ight  have  been  mot,  the  expenditure  in  scarce  resources  was 
greater  than  necessary.  As  systems  became  more  complex  and 
the  costs  associated  with  supporting  such  systems  astronom¬ 
ical,  it  was  evident  that  resource  expenditure  must  be  a 
variable  included  in  an  allocation  nodol.  Due  to  the  impor¬ 
tance  of  the  cost  consideration,  more  and  more  emphasis  has 
been  riven  to  this  quantity  in  allocation  prooeoduros.  This 
increased  emphasis  has  culminated  in  the  use  of  optiiuistion 
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techniques  for  seeking  and  defining  tho  most  rosource  effec¬ 
tive  methods  of  designing  and  developing  a  pioco  of  eauip- 
nent.  By  aoplying  optimisation,  a  best  alternative  can 
often  be  found  that  will  save  time,  effort  and  nonoy  in  the 
long  run#  A  brief  description  follows  of  the  more  popular 
optimization  methods  now  in  use. 

METHODS  OP  OPTIMIZATION 

Of  the  several  techniques  used  for  tho  optimization  of 
reliability  constrained  by  scarce  resources  the  more  popular 
methods  arc  Lagrange  multipliers,  linear  nrograjjjminn  and 
dynamic  programming.  As  would  be  expected,  each  of  those 
methods  varies  in  difficulty  of  formulation  and  solution 
according  to  tho  function  to  be  optimized  and  tho  nature  of 
the  constraint  relationships.  Here  an  attonpt  rill  be  made 
to  give  a  brief  explanation  of  how  each  method  is  employed. 

The  Lagrange  multiplier  method  is  analytical  in  that  it 
attempts  to  simultaneously  zero  all  the  partial  derivatives 
of  the  Lagrangian  function.  A  linear  combination  of  the 
objective  and  constraint  equations  is  first  formed,  then  the 
partial  derivative  is  successively  takon  with  respect  to  the 
variables  of  the  objective  function  and  the  Lagrangian  mul¬ 
tipliers.  The  resulting  simultaneous  equations  are  solved 
in  terns  of  tho  variables  of  interest  which  rill  thon  opti¬ 
mize  the  objective  function.  Two  principle  difficulties  are 
encountered  in  this  method.  The  first  is  the  difficulty 
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associated  with  the  differentiability  of  the  objoctive  and 
constraint  equations,  while  the  second  is  the  simultaneous 
solution  of  the  thus  obtained  equations. 

The  linear  programming  method  considers  objective  and 
constraint  functions  which  are  linear  in  nature.  The  liriear 
constraints  expressed  as  inequalities  are  made  equations  by 
the  use  of  slaelc  variables  and  are  arranged  in  matrix  nota¬ 
tion.  This  matrix  is  then  iterated  successively  (by  an 
operational  algorithm  suoh  as  the  Simplex  algorithm)  in  the 
direction  which  tends  to  optimize  the  objoctive  function. 
Other  solution  techniques  are  available  in  linear  program¬ 
ming  such  as  graphical  solutions  and  the  inverse  matr’x 
technique.  In  the  graphical  solution,  for  instance,  the 
region  of  feasible  solutions  (those  which  satisfy  t'~e  re* 
source  constraints!  is  first  found,  then  the  optimum  solu¬ 
tionis)  for  the  objective  function  is  examined  in  this  re¬ 
gion.  There  are  several  drawbacks  to  the  linear  programming 
model.  First,  the  objective  and  constraint  relationships 
must  be  linear  combinations  of  the  problem  variables  unless 
a  suitable  nonlinear  to  linear  transformation  is  found. 
Secondly,  the  easily  understood  graphical  tec'  niquo  is 
Ti~itod  to  three  variables  (three  dimensional  representa¬ 
tion).  Lastly,  matrix  operations  in  more  than  throe  or  four 
variables  are  quite  cumbersome  for  longhand  computation  and 
usually  require  digital  computer  solution. 

In  dynamic  programming  a  return  function  is  defined  as 


it  t.  , 
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the  sum  of  the  values  of  the  alternative  functions  of  return. 
The  value  of  each  of  those  alternative  functions  is  deter¬ 
mined  by  the  amount  of  resource  which  is  expended  on  that 
particular  task.  In  addition  it  is  permissible  to  have 
several  stares  at  which  portions  of  the  total  resource  can 
be  expended.  Thus,  tho  total  return  is  the  sum  of  the  return 
functions  at  each  stage  whore  resources  aro  expended.  The 
problem  is  then  to  optimise  the  total  return  function 
subject  to  the  resources  available,  iimploying  oxisvonce  and 
uniqueness  theroms  which  have  been  developed  it  is  possible 
to  make  successive  approximations  to  the  optimum  return  and 
approach  a  desired  optimum  policy.  The  main  problem  with 
this  technique  is  that  a  digital  computer  is  required  in  all 
but  the  simplest  cases.  Also,  the  number  of  numeric  itera¬ 
tions  required  in  problems  where  an  iterative  algorithm  is 
used  is  largely  dependent  on  the  closeness  of  the  initial 
approximations  to  the  optimum  policy. 

In  the  chapter  that  follows,  a  portion  of  the  litera¬ 
ture  available  on  solution  of  the  reliability  apportionment 
problem  is  reviewed.  The  work  in  these  several  papers 
considers  the  use  of  tho  aforomontionod  techniques  for  var¬ 
ious  cases  of  optimum  apportionment.  Chapter  III  considers 
the  optimal  allocation  for  a  series  system  where  the  system 
reliability  or  total  l:>  specified.  Two  cases  aro 

considered;  elements  with  a  constant  failuro  rate  and  ele¬ 
ments  v:ho36  failure  rates  oan  be  deseribed  by  a  Weibull  dis¬ 
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CHAPTER  II 


LITERATURE  REVIEW 

As  was  mentioned  in  the  introduction,  many  allocation 
methods  have  evolved  from  work  done  in  reliability  over  the 
past  two  decades#  The  scope  of  these  methods  ranges  from 
simple  to  quite  complex  depending  on  assumptions  made  and 
the  variables  included  in  the  development#  With  respect  to 
resources  expended  in  an  allocation,  the  mothods  so  far 
developed  can  be  divided  into  two  distinct  categories;  basic 
allocations  and  optimal  allocations# 

BASIC  ALLOCATIONS 

Basic  allocations  arc  so  called  because  the  model  used 
to  allocate  element  reliabilities  or  failure  rates  only  con¬ 
siders  the  overall  system  specification  to  be  met;  no  con¬ 
sideration  is  given  to  the  resources  exncnd^d  in  mooting 
this  snecificatios.# 

The  most  elementary  method  used  is  the  equal  allocation 
method  (1 )Y  Assumptions  for  this  method  are: 

1 •  independently  failing  components  and 

2*  serially  arranged  components. 

If  a  .:ystcm  reliability  goal  is  defined  as  R*  and  this  goal 

*  Number  in  parentheses  refers  to  numbered  references  in 
the  List  of  References. 
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has  not  boon  met,  then 

n 

R»  >  II  Rj 
i*1  1 

and  the  component  reliabilities  must  be  set  at 

Ri  “  ~y~*r~  • 

Although  this  method  is  straightforward  and  easily  applied, 
it  is  only  economically  suitable  when  tho  elements  are 
similar  enough  that  the  cost  of  reliability  increase  among 
them  is  nearly  equal.  If  this  is  not  the  case,  large  sums 
of  money  could  be  spent  in  increasing  the  reliability  of  a 
very  complex  component  when  it  could  be  spent  more  effec¬ 
tively  on  the  lower  cost  components. 

ARI1IC  (2),  has  developed  a  slightly  improved  allocation 
technique  by  defining  a  weighting  factor.  Assumptions  in 
the  development  of  this  technique  are: 

1#  independently  failing,  serially  arranged  components, 

2.  constant  component  failure  rates  and 

3.  initial  failure  rate  X  ,  is  known  for  each  com- 
ponent. 

The  weighting  factor  is  defined  as 


for  n  components  so  that  if  a  system  failure  rate  of  X  is 


desirod,  the  failure  rate  of  the  component  is 


Thus,  the  decrement  ^or  increments  of  eaoh  component  failure 
rate  is  proportional  to  its  initial  value.  Tuo  specific 
problems  arise  from  using  the  ARINC  method.  First,  if  a 
component  failure  rate  is  unreasonably  high  or  low  initially, 
the  allocated  failure  rate  will  not  tend  to  remedy  the  situ¬ 
ation.  Secondly,  no  consideration  is  made  for  the  effort 
to  decrease  a  components  failure  rate.  Less  total  effort 
could  be  involved  in  arriving  at  the  same  system  failure 
rate  if  other  component  failure  rates  wore  allocated. 

A  method  slightly  more  complex  is  the  AGR213  (7)  allo¬ 
cation  technique.  The  assumptions  for  tho  use  of  this 
method  are: 

1.  constant  element  failure  rates, 

2.  serially  arranged  components  where  the  itil  component 
contains  n^  elements  which  are  required  for  time  t^ 
in  each  mission,  and 

3.  the  ith  element  has  conditional  probability  that 
its  failure  x?ill  result  in  system  failure. 

If  overall  system  reliability  R*  is  desired,  the  mean  life 
allocated  to  the  oonponent  is 


n^-ln  d') 


n 

11  =  2 

i=i 


ni 


. 


where 
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This  method  is  considerably  more  realistic  than  tho  afore¬ 
mentioned  techniques *  It  not  only  considor3  tho  complexity 
of  the  component,  but  also  the  essentiality  and  mission  op¬ 
erating  time#  In  this  sense  it* does,  to  a  decree,  consider 
the  difficulty  involved  in  reducing  a  components  failure 
rate.  The  failure  rate  is  in  inverse  proportion  to  the  com¬ 
ponent  complexity.  If,  however,  the  elements  in  one  com¬ 
ponent  are  not  similar  to  the  elements  in  anothor  component, 
the  complexity  is  not  a  good  comparative  measure  of  tho 
difficulty  to  decrease  failure  rate.  Additionally,  it  is 
hard  to  dofine  a  components  essentiality  in  the  early  de¬ 
sign  phase. 

For  considering  redundancy  within  subsystems,  the 
1TAVW2P3  (6)  allocation  is  applicable.  Assumption  for  its 
use  are: 


1 .  subsystems  are  serially  arranged,  each  with  constant 
failure  rate,  and 

2.  the  failure  probability  of  each  subsystem  is  inde¬ 
pendent  of  other  subsystems. 

The  first  step  i3  to  reduce  the  redundant  configuration s) 
to  a  series  system  (i.e.,  any  active  element  results  in 
system  failure).  The  redundant  failure  rate  i3  estimated 


by 


\(t) 

SR  3 


yt) 


X 


s 


where  R  (t)  is  the  serial  pystem  reliability  over  timo  t, 

S 
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R  (t)  is  the  redundant  system  reliability  and  X„  is  the 
serial  system  failure  rate.  If  a  system  reliability  of 
R!(t)  is  required  the  reliability  of  the  i  subsystem  is 

Siv.n  by  Xi/X0 

Bjlt)  »  li'(t) 

where  X  ^  is  either  a  redundant  or  series  failure  rate  of 
a  redundant  or  series  subsystem  and  XQ  is  the  sum  of  all 
subsystem  failure  rates*  This  is  the  same  result  as  ob¬ 
tained  by  the  ARTEC  (2)  technique. 

As  with  the  other  basic  nethods,  no  consideration  is 
given  to  the  effort  to  incroase  a  susbsystem  reliability* 

In  addition  it  is  necessary  to  have  at  least  a  good  estimate 
of  the  series  and  redundant  reliabilities  for  a  particular 
time  as  well  as  series  failure  rate.  The  results,  as  with 
some  earlier  techniques,  are  directly  dopondent  on  prior 
knowledge  of  components  reliabilities* 

OPTIMAL  ALLOCATIONS 

Up  to  this  point  only  basic  allocation  techniques  have 
been  reviewed.  Nethods  which  will  now  be  discussed  are  not 
basic  allocations  in  that  parameters  eonsicterod  in  the  allo¬ 
cation  are  not  just  those  that  quantitatively  describe  the 
reliability,  ossontiality  or  complexity  of  the  systems  ele¬ 
ments.  Procedures  have  boon  developed  to  not  only  allocate 
in  so.:e  -artner  to  satisfy  a  reliability  requirement,  but  to 
do  t:iis  while  satisfying  some  constraint  on  rer,uuroo  expen- 
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dituros.  Several  of  these  methods  will  now  be  discussed. 

One  of  the  earliest  attempts  at  finding  an  optimum  al¬ 
location  with  a  minimum  of  resource  expenditures  was  devel¬ 
oped  by  Albert  (1 ) .  An  effort  function  is  assumed  that 
reflects  the  expenditure  (of  effort,  time,  money,  etc.) 
necessary  to  increase  the  reliability  of  a  component.  The 
assumptions  in  the  development  are: 

1.  serially  arranged,  independently  failing  components, 

2.  initial  component  reliabilities  are  well  established, 

3.  all  components  share  the  same  effort  function,  and 

l|>.  the  effort  increases  monotonically  with  component 

reliability;  an  increment  in  component  reliability 
requires  as  much  or  more  effort  at  higher  initial 
reliability  levels. 

With  these  restrictions,  the  minimum  effort  allocation  for 
an  n  component  systom  is 

1/k 


R 


R« 


m 

rx 

j«i 


where 


R=  reliability  allocated  to  components  with  initial 
reliability  less  than  'fTnT. 

R,js  systom  reliability  goal  • 

m 

**  reliability  produot  over  all  elements  with  ini- 

J*1  J  _ _ it 

tial  reliability  greater  than  m\Jr  (IT  R. 

v  j«1  * 

®1  when  no  components  meet  this  requirement) . 
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k=  number  of  components  that  have  reliability  less  than 

Where  a  system  is  encountered  that  meets  the  assumptions  of 
this  model,  the  allocation  procedure  is  simple  and  straight* 
forward.  The  assumption  of  a  common  effort  function  is 
rather  restrictive  in  that  this  will  probably  ocour  only 
when  the  components  arc  identical,  which  is  seldom  the 
case. 

For  considering  redundant  configurations.  Brier  (3)  has 
developed  an  optimisation  model  for  switched  components  in 
parallel  standby.  The  author  considers  a  system  with  N 
identical  subsystems  connected  in  a  parallel  switching 
arrangement,  ^ach  subsystem  is  oonposod  of  H  identical  ele¬ 
ments.  The  assumptions  made  in  tho  development  are 

1 .  all  elements  are  identical  and  have  constant,  inde¬ 
pendent  failure  rates, 

2.  elements  do  not  fail  when  not  energised, 

3.  successful  operation  of  one  element  is  adequate  for 
subsystem  operation, 

ii..  switches  between  subsystems  are  identical  with  time 
independent  probability  of  suooess,  and 

5.  switch  failure  is  duo  only  to  non-switching  when 
inquired;  premature  switching  is  not  considered. 

The  s  vs  torn  begins  operation  on  the  first  subsystem  until  all 
!:  elements  have  failed  then  switches  sequentially  until  all 
1!-1  redundancies  are  failed  or  a  switch  fails.  An  optimum 
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N  and  K  are  found  from  the  Lagrangian 

If  H  - 

L  =  (1-?)  2  T  -  o-H  (1+rM) 

i-11 

whore  P  is  tho  switch  reliability,  cr  is  a  Lagrange  multi¬ 
plier,  and  r  is  the  ratio  of  element  cost  to  switch  cost* 
Prom  algebraic  rearrangement  the  system  mean  time  between 
failure  can  be  found  from 

1  N-1  &  1 

I7?BPayst.em  -  ^  (1-P)  9 

where  1/X  is  the  HTBP  of  each  element.  The  total  cost  is 
found  from  the  second  part  of  tho  Lagrangian  function  and  is 


where  C  is  the  cost  of  a  svritch.  Sinco  a  closed  form  solu¬ 
tion  was  not  obtainable  for  the  above  Lagrangian,  the  author 
provides  a  digital  computer  program  and  shows  a  table  of 
optimum  redundancies  for  ranges  on  the  problem  variables  of 
1  ft  i;  ft  50,  1  ft  K  *  50,  .95  -  P  -  .99  and  .25  *  r  *  2.00. 

The  cost  model  developed  is  a  poworful  tool  to  the  designer 
considering  a  switched  parallel  redundancy#  Several  draw¬ 
backs  exist  in  this  model: 

1 •  Failure  rntes  must  bo  constant. 

2*  All  elements  and  switches  must  be  identical. 

3.  It  is  necessary  to  have  access  to  the  authors  com¬ 
puter  program  or  tables  of  optimum  redundancies. 

A  field  very  much  akin  to  reliability  is  availability# 
Availability  considers  not  only  reliability  or  failure  rate 
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of  an  equipment,  bnt  also  the  time  required  to  brine  t'ie 
equipment  from  a  failed  state  to  an  operable  state# 

IlclCichols  and  Tessor  (  )  have  *onsiderod  the  dn •»«*/.*  tlon  of 
cost  to  achieve  a  level  of  inherent  availability*  Although 
not  oxprossly  designed  for  use  as  a  reliability  optimization 
technique,  a  form  of  the  equation  is  presented  that  can  be 
used  as  a  tool  for  failure  rate  allocation#  The  assumptions 
of  the  model  are 

1  *  serially  arranged  components  with  independent  fail¬ 
ures, 

2.  each  component  has  a  constant  failure  rate,  and 

3*  the  cost  of  a  decrement  ax  in  clement  failure  rate 

>1 


is  a 


I 


Inherent  availability  can  be  expressed  as 


where  is  a  constant. 


A  = 


i  +  2  Mj  X  j 

i=i  1  1 


for  n  elements  in  series  where  is  the  mean  time  to  repair 
the  itn  item,  and  is  the  failure  rate  of  the  Item# 

If  wc  let  the  equal  unity  for  oaoh  element,  the  inherent 
availability  reduces  to 


n 

i  +  2  Xi 

1=1  1 


n 


•or  a  sot*io3  avatar.,  however,  j  \  i5  the  system  failuro 

1=1  1 

rate#  Usinq  a  Laqranqian  of  the  form 
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where  A  is  the  present  availability,  A1  is  the  desired 
availability  (determined  in  this  case  by  1/(1+ X*),  A*  being 
the  desired  system  failure  rate),  and  K  is  the  Lagrangian 
multiplier*  ttearrangini?  the  results  of  the  author,  the  min- 
inun  cost  failure  rate  allocation  for  the  i*  element  is 


This  allocation  is  considerably  more  flexible  than  that  de¬ 
vised  by  Albert  (1)  in  that  each  element  can  have  a  differ¬ 
ent  cost  function.  Two  possible  drawbacks  arc  the  necessity 
of  a  series  as  sumo  ti  on  and  the  aocuracy  of  “fitting  1  the 
cost  function  given  to  the  actual  cost  versus  failure  rate 
data.  The  cost  function  proposed  for  the  scries  system  in 
the  next  chapter  is  more  powerful  in  that  it  is  a  two  para¬ 
meter  family  of  curves  and  can  therefore  nore  accurately 
describe  a  general  cost  versus  failure  rate  curve. 

A  paper  written  by  3pivey  (0)  compares  and  contrasts 
the  results  obtained  by  using  three  different  optimization 
methods  on  a  specific  system  configuration.  The  nothods 
considered  arc  dyna  mic  prorramming,  mathematical  pro~rnrring 
arid  Lagrange  multipliers.  Due  to  the  complexity  of  the 
system  considered,  a  series-parallel  notworlc,  the  author 


elected  to  employ  a  digital  computer  solution.  The  spocific 
model  used  with  each  technique  will  not  be  repeated  here, 
only  the  results  found.  It  was  determined  that  tho  optimum 
solution  for  the  system  considered  was  found  by  dynamic  pro¬ 
gramming.  The  solutions  obtained  by  mathmatical  programming 
and  Lagrangian  multipliers  were  within  ,6rf  and  ,5: '  of  the 
optimum.  These  errors  night  be  attributiblo  to  arithmetic 
operations  in  the  digital  computer  or  tho  nature  of  the 
model  developed. 

The  following  chanter  will  discuss  and  develop  a  sories 
model  and  some  constraints  applicable  to  tho  allocation 
problem.  Por  the  series  configuration,  two  forms  of  the 
cost  function  are  considered.  The  preference  of  one  devel¬ 
opment  over  the  other  should  be  based  on  tho  adequacy  of  a 
particular  cost  function  in  describing  actual  improvement  of 
reliability  or  failure  rate  with  cost.  In  Chapter  IV  an 
optimum  allocation  is  developed  for  the  parallel  system. 


CHAPTER  III 


SERIES  SYSTEMS 
INTRODUCTION 

In  this  chapter  only  series  systems  will  be  considered. 
This  configuration  is  a  result  of  the  following  assumptions: 
1 •  The  probability  of  failure  of  any  element  is  inde¬ 
pendent  of  failures  in  any  other  elements. 

2 •  All  elements  are  required  to  operate  for  system 
success. 

The  first  two  sections  deal  with  a  series  system  with  one 
additional  assum »tion,  constant  element  failure  rates.  A 
specific  cost  function  is  developed  for  which  the  optimum 
allocation  is  found  for  both  a  total  cost  constraint  and  a 
system  failure  rate  specification.  Constant  element  fail¬ 
ure  rates  are  not  assumed  for  the  second  series  model  con¬ 
sidered. 


co:;sta::t  failure  rates  -  total  cost  constraint 

A  special  case  of  the  series  model  is  where  oach  ele¬ 
ment  displays  an  exponentially  distributed  time  to  failure. 
V.hon  this  is  true  and  each  element  failure  is  independent 
of  all  others,  the  statement 


is  ecui valent  to 


m 

maximize  11  R. 

i=1  1 

m 

minimize  £ 
i=1 


where  is  the  failure  rate  of  the  i  element  and  is 
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tho  reliability  of  the  i^  element.  In  this  model  the  ob¬ 
ject  will  be  to  minimize  the  system  failure  rato  subject  to 
an  overall  expenditure  in  achieving  this  failure  rate. 

Thus,  it  is  necessary  to  express  tho  cost  of  each  oloment  as 
a  function  of  the  failure  rate  achieved  in  that  elenent. 
Sunning  over  the  costs  of  these  elements,  tho  total  system 
cost  can  bo  determined. 

One  reasonable  assumption  in  expressing  the  cost  of  an 
element  as  a  function  of  its  failure  rato  is  that  addition¬ 
al  expenditures  for  failure  rate  improvement  result  in 
diminishing  returns  of  failure  rate  reduction.  One  func¬ 
tion  that  displays  this  characteristic  is  an  exponential 
cost  function  in  which  the  cost  increases  exponentially 
with  a  reduction  in  failure  rate.  Expressed  mathematically, 
the  cost  of  failure  rate  in  elomont  i  is 
Bi(e0l-ei) 

A^  e 

where  is  the  cost  of  obtaining  the  lowest  possible  fail¬ 
ure  rate  eo^  in  element  i,  and  is  tho  cost  gradient.  To 
completely  define  this  cost  function,  two  estimates  are 
required;  1 )  the  cost  A  at  tho  minimum  failure  rate  obtain¬ 
able  and  2)  a  cost  £  at  the  present  or  any)  failure  rate 
level  For  a  cost  C  at  failure  rate  level  9,  B  can  be 
estimated  as 

3  =  ln(c/A)/(eo-S) 

With  a  cost  function  defined  for  each  element  of  the  series 
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system,  the  total  cost  for  m  elements  is 
i«1  1 

If  for  the  development  of  this  system  the  upper  bound  on  ex¬ 
penditure  is  H,  the  Lagrangian  function  can  be  written 

8.  via  M®«<-®<)  1 

L  «  2  -  X  $  As  -  1  ol  1  "I 

i=1  x  [i=1  ^ 


Talcing  partial  derivatives,  the  equations  to  satisfy  are 
1  +XAj_B^  e  =  0  i=l,2,...,m 


£U 

and  £  A.  a 


Bi<®oi-V 


-  H  =  0. 


One  additional  restriction  is  that 


9i  *  9oi 


isi  i^i  « jin 


to  insure  that  a  failure  rate  is  not  allocated  that  is  bet¬ 
ter  than  can  be  obtained*  Sunroing  Equations  1  over  all  m 
olenents  gives 

I  J  .xl  h  .  o. 

i»1  Bi  i=1  1 
Substituting  Equation  2, 

\=  1  8.  1 

-  5  2  sr  • 

H  i-1  ®i 


Kq,  1 


£q.  2 


l'sin~  this  in  aquation  1  and  solving  for  9^,  the  optimum 
failure  rate  allocation  for  the  1th  element  is 
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®i  =  ®oi  +  \ln 


in 

AiBi  X 
1  1  1=1 


In  the  case  Where  this  last  inequality  is  not  mot  by  tho 


allocation  for  one  or  more  elements,  the  failure  rato  of 


these  elements  is  set  equal  to  their  respective  minimum 


failure  rate.  This  will  occur  when  the  quantity 


Aj3.  X 
1  1  1=1 


is  found  to  be  less  than  one  in  a  particular  elemont.  At 


this  point  a  reallocation  must  be  mnde  to  the  remaining  ole-* 


nents  with  the  remaining  resources  available  for  allocation* 


The  new  maximum  expenditure  is 


Hnew  =  H  -  XAt 


for  the  remaining  elements,  where  the  summation  is  over  the 


k  elements  sot  at  their  minimum  failure  rate. 


iSXAHPL3  1 


During  the  design  phase  of  a  vibration  recording  sys- 


of  failure  rate  to  the  system.  The  system  is  composed  of 


accelerometer,  amplifier  and  recorder.  All  elements  must 


function  for  3ystom  operation  (series  assumption),  and  ele¬ 


ment  failures  are  independent.  The  elements  are  assumed  to 


have  a  constant  failure  rato  and  the  applicable  cost  para¬ 


meters  are: 


Sq.  3 


ten  it  is  desired  to  find  the  most  cost  effective  allocation 
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Accelerometer(1 ) 

Amplifier (2) 

Recorder (3) 

®o 

.0001 

.001 

.002 

A 

•<>250 

$1600 

$1400 

'e 

.00025 

.003 

.007 

3  ! 

?ioo 

$900 

.3500 

The  total  expenditure  for  system  development  should  not  ex¬ 
ceed  $2,500.  What  failure  rates  should  be  allocated  to  the 
system  elements? 

Using  ?  and  *c  for  each  element,  the  cost  gradient  3^ 
can  be  found  as  1^=6110,  B2=346#  02=199#  Tho  initial  allo¬ 
cation  using  Equation  3  is  8^=  .000365,  .00308  and 

6y  .00146.  At  this  point  the  system  failuro  rate  is 
. 004905.  Howevor, .  since  6^  less  tilan  obtainable,  it 
is  set  to  its  lowest  value  of  .002.  This  results  in  an  ex¬ 
penditure  of  .,>1,400  on  element  three.  Using  the  allocation 
equation  again  with  Hnew~  $2,500-  ul  ,400  «  $1,100  gives 

8^  .000337  and  e2=  .00258  . 

The  design  criteria  for  minimum  system  failure  rate  is 


Accelerometer  Amplifier  Recorder 


e 

.000337 

.00258 

.002 

Cost 

•958.60 

?i  ,041 .40 

31 ,400 

System  failure  rate 

=  .004917 

C0ITSTACT  FAILURE  RATES  -  SYSTEM  FAILURE  RATE  CO”;?  tAIOT 


Tho  problem  considered  here  is  essentially  tho  some  as 

\ 

I 

. — — nwnfrjTt  nnr  jf  m 


the  one  just  presented  except  that  in  this  c aso  allocation 
must  satisfy  some  requirement  on  system  failure  rate.  Stat¬ 


ed  mathematically. 


minimize  2  A,e  x  01  x 
i=1  1 


m 


subject  to  £  6,  *  e.  and  6,  »  e_  i=1 ,2 
i=1  -  f  is 

The  Lagrangian  function  can  be  written  as 


»•*! 


m. 


m  3.  (8  *  ) 

L  =  2  A,e  1  01  1 
i=1  1 


-  X(  2  e.  -  e  ) 

i=1  1  8 


where  the  conditions  to  satisfy  in  this  caso  are 
3,(e  .-e,) 


+\  = 


i^i §2,  •  •  •  • , m 


and 


2  - 

i=i  1 


6S  -  0. 


Taking  the  natural  logarithm  of  Equation  k  and  summing  over 
all  i  yields 

*L  S  \  »  1  n  1 

2  *  2  6oi  -  ln(-X)2  5  +  2  ?  InfA^B,)* 

i=i  1  i=i  01  i=iBi  i=iBi  x  1 

Substituting  iquatlon  J>  into  this  equation  and  solving  in 

terms  of  the  Lagrange  multiplier  give 

rf'  m 


ln( 


"i,  VS»  1 

.  2  e  +  2  i  in(A.B, )  -  e 

-A)  =  igl  01  i=13i  1  1 _  s  , 


m  , 

2  g 

i=iBi 


Using  Equation  Ij.  with  this  relationship  gives  the  optimal 
failure  rate  allocation  to  the  i  element  as 

is1 , 2 , • • ,m 


®i  -  ®oi  +  lL  [ln<Ai3i)  - 


2qs. 


5 


j£q.  6 
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and  0i  ^  8q^  1=1, 2,..., m  • 

Again,  if  this  laat  inequality  is  not  net  by  the  first  allo¬ 
cation,  then  for  those  elements  for  which  it  i3  not  satis¬ 
fied  the  minimum  failure  rato  should  be  allocatod.  ?hi3 
situation  will  occur  whenever  the  quantity  A.3^  is  less  than 
-  X  #  With  these  elements  allocated,  the  new  system  failure 
rate  to  satisfy  is 


8, 


new 


=  e 


summed  over  all  k  elements  at  their  minimum  failure  rate. 


EXAKPLfci  2 

For  the  system  presented  in  Sample  it  is  decided 
after  further  investigation  that  a  system  failure  rate  of 
.00$  would  be  adequate  to  moot  design  specifications,  What 
amount  should  bo  spent  on  each  element  for  development  and 
what  failure  rates  can  bo  expected? 

Using  the  B^  and  8oi  as  defined  in  Example  1, 
ln(- A  )  =  12  .61;.  This  is  greater  than  ln(A0B^)  so  that  8^ 
is  again  set  to  .002.  ‘±ho  now  failure  rate  constraint  ie 

58  .00$  -  .002  =  .003*  Allocating  .003  between  the  re- 
naininr;  elements  yields 


61  =  .000330  and 

'2ho  new  design  criteria  are 


8 >  =  .00267. 


Accelerometer 

Anplifiop 

Recorder 

®1 

.000330 

.00267 

.002 

Cost 

361.30 

?1,010 

>i  ,400 

-  *^****t^Jai*asi*  i+f  ■  ■ ffiaaic. &ta* *  ih  f  -aeAwadi  r <£. 


2k 

The  system  failure  rate  is  now  .005  at  a  total  cost  of 
&,! 4.71.30. 

Because  many  systems  exist  which  do  not  exhibit  a  con¬ 
stant  failure  rato  the  section  which  follows  will  consider 
a  series  system  where  the  elomonts  can  have  othor  than  a 
constant  failure  rate.  It  is  assumed  that  each  olonents 
orobability  density  function  can  be  adequately  described  by 
the  Ueibull  distribution. 

A  H03B  GENERAL  SERIES  MODEL  —  DISCUSSION 

Up  to  this  point,  the  only  series  system  considered  was 
one  which  consist ed  of  elomonts  ’which  all  displayed  exponen¬ 
tially  distributed  times  to  failure.  The  optimum  allocation 
derived  is  only  valid  where  the  failure  rate  for  each  ele¬ 
ment  is  constant  and  the  form  of  the  cost  function  assumed 
can  adequately  describe  the  cost  of  achoiving  a  failure  rate 
for  an  element.  Even  though  the  constant  failure  rato  as¬ 
sumption  i3  not  too  far-fetched  for  most  electronic  equip¬ 
ments,  it  is  definitely  not  applicable  to  contain  other 
equipments.  In  practice  it  is  found  that  certain  equipment 
nay  display  tines  to  failure  that  are  distributed  as  a  nor¬ 
mal,  gamma,  chi-square  or  othor  familiar  probability  density 
function  (p.d.f.).  For  such  equipment,  it  is  believed  that 
money  and  effort  spent  in  design  and  development  will  result 
in  a  more  reliable  piece  of  equipment.  This  being  tho  case, 
it  should  be  possible  in  certain  instances  to  find  tho  basic 
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relationship  between  the  levol  of  reliability  achieved  and 
the  effort  expended  in  the  process.  In  the  case  of  the 
exponential  p.d.f.  it  is  possible  to  relate  the  reliability 
to  cost  by  defining  a  functional  relationship  botwocn  the 
failure  rato  and  its  associated  cost.  Such  a  function  was 
defined  in  the  allocation  method  dorivod  in  the  first  part 
of  this  chapter.  Generally,  for  a  p.d.f.  other  than  the 
exponential  or  chi-square,  it  is  nocossary  to  specify  more 
than  on©  parameter  in  order  to  specify  the  distribution. 

The  resulting  functional  relationship  between  cost  and  re¬ 
liability  in  an  element  t*ill  bo  more  complicated  since  it 
may  now  be  a  function  of  all  the  distribution  parameters. 

At  this  point  if  the  discussion  is  restricted  to  the  fam¬ 
ily  of  two  parameter  p.d.f^s,  a  co3t  function  for  obtaining 
distribution  parameters  Cl  and  (2  would  be  0(0.,/?).  An 
example  of  such  parameters  might  be  the  mean  and  standard 
deviation  of  the  normal  distribution. 

Under  the  assumptions  of  a  series  system,  the  system 
reliability  is  the  product  of  the  element  reliabilities. 

The  system  reliability  is  thus 

lUt)  =  TI  Rj(t) 

3  i=1  1 


for  n  elements  and  t  hours  of  operation.  This  is  equivalent 
to  the  system  failure  rate  being  the  sum  of  the  element 
failure  rates,  or 

i 

i 

A 


•tf 
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XL 


^  (t)  «  2  ®*(t)  # 
3  i*i  x 


For  any  p.d.f.,  the  reliability  can  be  written 

\(t)  =  exp  -J%  e^ujdu 

th 

where  9^{t)  is  the  failure  rate  of  the  i  element*  In  or¬ 
der  to  maximize  the  system  reliability  at  time  t  it  Is  ne¬ 
cessary  to  find  the  minimum  value  of 


Jo  'Jt  [IaH  * 


Sq.  7 


3q.  8 


Sq.  9 


for  the  system  elements.  To  greatly  simplify  the  task  of 
minii:iisinR  Equation  9  with  respect  to  tho  distribution  par¬ 
ameters  of  each  element,  a  common  and  intolerable  p.d.f. 
could  be  used  for  each  element.  A  p.d.f.  which  is  quite 
flexible  in  its  ability  to  approximate  othor  distributions 
is  tho  V/eibull  distribution  which  can  be  expressed  as 

n  n.  Z?-1  -v/a 
f(t;a,/?)  =  (/i/a)  t  e  ,  tso. 

For  judicious  choices  of  a  and/?,  acourate  approximations 
to  tho  normal,  exponential,  gamma,  chi-square,  and  other 
distributions  can  be  obtained.  Ihe  Appendix  considers  the 
values  of  CL  and  that  can  bo  estimated  to  describe  other 
p.d.f.’s.  Since  the  failure  rato  of  the  V/eibull  distribution 
can  be  expressed  as 

•1 


9(t) 


R  (S- 1 
fr  * 


io.  10 


the  right  hand  side  of  Equation  9  can  be  written 


1 1  lit 
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j-:  iL-H « <  [|,«  ■ft> 

which  when  evaluated  yields  for  the  system 

r*  rs  .tWi .  t  £ 

J  0  L  i=1  J  i=1  ai 

Now  that  the  integral  of  the  system  failure  rate  in  the 
interval  of  interest  (0,t)  has  been  expressed  as  a  function 
of  the  tine  interval  and  Weibull  parameters  for  each  ele¬ 
ment,  it  i3  of  interest  to  define  the  costs  associated  with 
designing  and  developing  the  VJeibull  parameters  CL  and  (2 
into  each  element.  Although  there  are  many  mathematical 
possibilities  for  expressing  tho  cost  of  achieving  parameters 
CL  and  Pin  an  element  as  a  fraction  of  these  parameters,  it 
might  be  helpful  to  understand  how  CL  and  P  affect  the  moan 
and  variance  of  time  to  failure  for  the  system.  It  can  be 
seen  fr on  Equation  11  that  each  elements  contribution  to 
tho  value  of  the  integral  Increases  with  increasing  P  and 
decreases  with  increasing  CL  .  Qualitatively  thon,  increas¬ 
ing  CL  and  decreasing  P  for  any  element  will  result  in  a 
smaller  value  of  the  integral  in  Equation  11  and  thus  a 
higher  system  reliability.  examination  of  the  equations  in 
the  Appendix  reveals  the  changes  in  the  distribution  mean 
and  variance  with  changes  in  CL  and/?.  It  is  intuitive  that 
for  a  given  element  the  cost  of  development  will  increase 
if  a  is  increased  and/or/?  is  decreased  in  tho  ./oibull  den¬ 
sity  describing  the  elements  p.d.f.  of  tine  to  failure. 


£q.  1 1 


► 


A  cost  relationship  which  is  proposed  for  use  in  the  follow- 
ing  analyses  is 

c(a  tj3)  =  Aa/?'b 

where  A,  a  and  b  are  positive  constants  which  appropriately 
describe  the  cost  of  a  particular  element.  Equation  12  is 
flexible  enough  that  a  v;ide  variety  of  actual  cost  versus 
reliability  curves  could  be  approximated.  The  Appendix  con' 
aiders  methods  to  approximate  parameters  in  thi3  co3t  func¬ 
tion  for  an  elemont.  One  such  sot  of  cost  curves  for  cost 
divided  by  A  versus  element  reliability  is  shown  in  Figure 
for  b=1.C  and  a=.25. 


whero  H  is 


ALLOCATION  OP  a  AMO  /3  POM  A  SfflllSS  SY5T-.M - 

TOTAL  COST  RiSqUIR-iLuTT 

Using  the  cost  function  proposed  in  the  last  section, 

the  allocation  problem  for  the  series  system  is  to  minimise 

n  +Pl  n  a,  --b, 

t  o  subject  to  %  A,  CL,X  *  =  II, 

i=i“o7  i=i  1  1  1 

the  maximum  permissible  cost  and  n  is  the  number  of  system 
elements.  The  Lagrangi&n  function  for  this  problem  is 
n  4,/^i  N  rn  a.  n 

L-h\i  -»]  • 

Taking  partial  derivatives,  the  equations  of  interest  are 

N  0-b  ± 

Pi  -  0  i=1,2,..,n 


a 
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ilq.  1 2 


F.qs  .13 


COMPONENT  PELIABtLITY 


4  V  “I 

ln(t0)  +  Xb^A^GL^  (3 ^  »  0  is1»2#««»n 


2qs*  14 


^-b4 


**  “4 

and  X  A,  CL/j&  1  -  H  *  0  • 
i=1  1  1  1 

Equations  13  and  M\.  oan  be  combined  to  give 


"  at  la(t0)  • 

Substitution  of  Equation  15  into  Equation  13  or  lit  yields 


ai  =  ih. 


a4 

or  O.^  =  -«» 


&iln(t0  j, 


Xa^A^ 


apT~ 


The  cost  constraint  can  nov/  be  written; 


H  B  4  -  *r  I  A 

i«1  l A/ 


•b^t  i  a£l  1 


ilq.  15 


Kq.  16 


VT  do.  17 


2quation  1 7  can  be  solved  for  the  Lagrange  multiplier  X  by 
trial  and  error  choice  of  the  multiplier*  Once  a  suitable 
multiplier  is  determined  it  can  be  used  in  Equation  16  to 

i.1. 

find  the  allocation  of  CL^  to  the  i  olement*  If  an  ele¬ 
ment  is  allocated  an  CL  and  (3  which  is  better  than  can  be 
obtained  for  that  element,  the  CL  is  set  at  the  highest  and 
(2  to  the  lowest  possible  values  for  the  element*  A  now 


allocation  is  now  made  over  the  remaining  elements  with  the 
renainin-  funds. 
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To  find  the  appropriate  multiplier  in  Squall  ••a  17  with  any 
deer©0  of  accuracy  would  require  a  digital  compufcor  solu¬ 
tion,  A  situation  where  the  solution  of  Equation  17  for  the 
multiplier  is  simplified  is  considered  in  the  following 
section. 


OPTII -SUM  ALLOCATION  P03  &7UAL 
hnen  a  situation  is  encountered  in  which  the  a^*s  can 
be  set  to  a  common  value  in  the  element  cost  functions, 
Equation  1?  can  be  solved  explicitly  for  the  multiplier. 
With  a^  «  a  i=1,2,...#n#  Equation  17  becomes 


a 


Equation  16  can  now  be  written 


Fran  Equation  Ifi  the  /X  can  bo  written 


Sor.o  mention  should  ho  made  of  the  situation  in  which  the 


3q*  1 6 


2q.  19 


a. *s  are  equal.  Ono  obvious  instance  is  where  the  p.d.f* 
of  tir.o  to  failure  for  all  elements  is  similar  enough  that 
the  ratio  of  cost  increase  to  increase  in  Cl  is  approximately 
equal.  It  is  possible  hotiover  for  ensos  to  occur  ’..hero  the 


32 


a£*s  are  nearly  equal  although  tho  p.d .f.  of  time  to  failure 
for  elements  is  not  t'^o  same.  The  example  which  follows 
will  illustrate  an  approach  to  the  case  w;erc  the  p.d.f.  of 
element  time  to  failure  is  not  the  sane  but  the  a^fs  aro 
nearly  equal. 


MXAKPLS  3 

It  is  desired  to  find  the  optimum  distribution  paramet¬ 
ers  to  allocate  three  elements  which  comprise  a  series  sys¬ 
tem  in  order  to  maximise  the  system  reliability  for  20  hour  * 
operation  with  the  funds  available.  Tho  p.d.f. fs  of  time  to 
failure  for  the  first,  second  and  third  elements  aro  normal, 
gamma  and  exponential.  Two  estimator  or  cost  versus  mean 


lifo  for  the  elements  are  as  follows: 

2  .2 


riornal 

<r  /  =  .01 

^  =  5 0  hrs. 

O 

II 

O 

O 

O 

2  2 

fi2  =  75  11 

C2  «->3  9  000 

Oanma 

<r  /  (i  =  .0*025 

* 20  " 

C1  =.;500 

0  0 

0 

rn 

» 

C\I 

C2  =.’jf.00 

aeponen.  &  i  fl/=  1*°  /i-1  ■  “>00  ;1  C1  =;!i200 

fl2  -  150  "  C2  =  ?220  . 

'..'hat  CL  and  /?  r^nmetarn  should  bo  allocatod  each  element  if 


total  cost  should  bo  )1,500? 

"  siny  the  methods  of  the  Appendix  to  estimate  tho  cost 
par  a:  rotors,  the  cost  functions  for  elements  one,  two  and 
throe  can  now  bo  writton 

16.47  .226  ,0-8.20  8.20  .228  ->-3.4‘: 


CL  *  *  Q  ~ 

1  » 


8.20  .228  ^-3.4‘: 

•  a2  Pi 


1 1  .nr 


■*****s® 
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and 


4.25  .220  O-.605 

e  3.,  fj. 


3  ^3 

Osin.'?  Equation  1C,  the  optimal  CL^  are  found  to  be 


a,  =  e 


40.55 


a2  »  o 


16.15 


CL  ^  s  0 


7-137 


.substituting  those  into  their  appropriate  C03t  equations, 
the  amount  to  be  srent  on  components  one,  too  and  three  is 


Ct  ~  :*>202  C2  a  39J4.O  C3  a  ^360  , 

for  a  total  cost  of  31, $02.  The  resultant  mean  lifo  of  each 
element  will  be 

Mi  =  31.2  hra.  fl2  -  34.7  hra.  =  1,250  hrs. 

The  system  reliability  can  bo  written  as 

|  tA 

exp  -X  -75—  f 
i=1  ai 


so  for  this  allocation  and  t  «  20  hours,  the  avsten  relia- 

o 

bility  is  #933*  An  improved  allocation  could  now  be  accomp¬ 
lished  by  making  now  cost  parameter  estimates  in  the  region 
of  the  initial  allocation. 

The  section  which  follows  oonsiders  a  problem  very  sim¬ 
ilar  to  that  just  discussed.  How,  however,  the  ays  tom  reli¬ 
ability  will  be  held  constant  and  minimum  003 t  found  to  sat¬ 
isfy  this  reliability. 


:  cg.tp  so*  svots:  *kmadility  ?  :  t  ?  j  - 

The  last  section  considered  maxird^iny  system  reliabil¬ 
ity  subject  to  a  total  expenditure.  A  typical  problem  anal¬ 
ogous  to  this  is  finding  the  minimum  design  and  development 


■j? 


2k 


cost  associated  with  mooting  the  system  reliability  require¬ 


ment.  Suppote  a  system  reliability  requirement  of  R  (t  )  is 

s  o 

specified  for  operating  time  tQ.  For  a  scrios  system  then 

n  .  Pi 


W  =  exp  -2  to 

i~i  a 


so  that 

a,  ,A 

S.-fc  ' ' ln  W  • 


In  the  following  discussion,  tho  tom  -  In  i\(tQ)  will  be 


considered  the  requirement  to  moot  and  will  be  denoted  by 


K#  Therefore  ^s(tQ)  is  equal  to  exp(-X).  ?or  tho  system 


composed  of  n  elements  in  series  the  Jagrangian  function  is 


2,  ai  n~h i  >rn  1 

L  -  2  a  a^p.  -X  2  1°  -  k 

i=i  1  1  1  Li=imrr  J 


where  the  object  is  now  to  minimize  the  total  cost  torn 


u  •  w“b. 


iw#:1 


The  equations  resulting  f^om  partial  differentiation  are 

ai-1  /0"bi  4.^ 

aiAiai  Pi  +  XI°  =0  i=1,2,..,n 


a,  _  -b .  -1 


Vi01!1  Pi  1  +  xi°  (t0)  i»1 n 


2  to  -  K  =  0. 
i=<  a. 


^qs*  19 


£qs.  20 


,  -  j-‘.y  i‘riiiA'T  rrf  lV'‘rAMrir~y)'i-*-:fc-*^  --  -'.-..Mb 
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Solution  of  bqualions  19  and  20  for  j3 i  gives 

n  _  bi 

-  r~TnTt0T 

Pron  this  it  is  seen  that  the  optimal  level  for  P in  an  ele¬ 
ment  is  the  same  whether  optimizing  in  terms  of  total  cost 
or  system  reliability  and  is  only  a  function  of  the  cost 
parameters  a  and  b  and  the  operating  time  tQ.  Equation  19 
can  be  used  with  the  equation  for  P^  to  find  the  CL.  as 
r.l  b<  |b4  b4/a<-|l-rv 


ai  Ai 


The  reliability  constraint  can  now  be  written 

I  ufn*  (.V‘1)^T  [MM  bl 

i='  L  ai  A±  J 


=  K. 


the  nultiplior  is  determined,  the  CL^  can  bo  found  fr^m 

,  n  bt  |  b4  bi/nq  1 

a,  =.  U)  v*  Iqisu7)  •  °i  ’ '  . 


L  ai  Ai  J 

The  following  section  doals  with  the  case  of  equal  a^  and 
the  resulting  straightforward  solution. 

ALLOCATION  FOX  3PUAL  ax 

Inspection  of  Equation  21  reveals  that  when  the  a^  are 
equal  the  CL  can  be  found  directly.  Letuing  all  a^  equal 


<2q.  21 


Prom  this  relationship  a  trial  and  error  teohnique  can  be 
used  to  find  the  required  multiplier.  As  before  an  accurate 
allocation  would  require  a  digital  computer  solution.  Once 


•*.q.  22 
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the  common  valuo  a,  Hquation  21  can  be  written 
-  -  -1  K 


n  b,  /a+1 

* 

2  e  1 
i=1 

[ailnrt0)j 

- 

a 

mJ 

aTT 


The  allocation  for  CL  and  J2  in  the  ith  ©lomont  becomes 


ai  - 


a“Jblt0) 


i  b>i  b^ai 
e 


a  Aj 


1 

ia+T 


b^/a+1 


i=1 


e 


bn-1 


|almt0)J 


a  A. 


a+T 


K 


A- 


a  ln(tQ) 


Sqs.  23 


aXAKPLH  ij. 

?or  the  system  of  Example  3  it  was  found  that  the  allo¬ 
cation  resulted  in  a  system  reliability  of  ,933  for  20  hours 
operation  at  a  cost  of  H ,  300.  ’.'/hat  savings  could  be  made 

if  the  system  reliability  requirement  for  20  hours  operation 
was  .90  (Ks.1052)? 


Using  the  cost  functions  as  developed  for  elements  one, 
two  and  three,  liquation  23  gives  the  following  Cl  allocations: 


ai  =  e 


WM3 


a2** 


17.73 


a  =  e 
3  6 


6.717 


These  valuos  yield  a  system  reliability  of  ,900,  The  re¬ 
sultin'*  element  oosts  are: 

c1  *  iidi  c2  =  $050  c3  =  ^325  • 


The  no w  total  cost  is  >1,356;  a  savings  of 

Analysis  of  the  type  just  presented  oould  be  a  valuable 
tool  to  examine  cost  effective  reliability  specifications • 

i 


KwwHMunmw 
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A  trade-off  could  be  performed  to  determine  a  suitable  bal¬ 
ance  between  cost  and  achieved  reliability*  In  this  way  a 
more  effective  system  could  be  defined. 

Chapter  3  has  restricted  the  discussion  to  scries  sys¬ 
tems,  The  first  part  considered  the  allocation  when  the 
elements  displayed  a  constant  failure  rate.  The  second  part 
of  this  chapter  considered  the  series  system  where  the  fail¬ 
ure  rate  of  each  element  was  not  necessarily  constant,  but 
the  distribution  of  times  to  failure  for  each  element  could 
be  adequately  described  by  the  Weibull  distribution.  An 
optimum  allocation  for  a  parallel  system  will  be  treated  in 
the  next  chapter.  The  VJeibull  approximation  will  o^ain  be 
used  to  aid  in  the  mathematical  development. 


CHAPTER  IV 


PARALLEL  SYST13MS 


INTRODUCTION 


This  chapter  will  bo  devoted  to  tho  consideration  of 


parallel  systems.  For  tho  purpose  of  discussion,  a  parallel 
system  will  be  defined  here  as  one  which  requires  the  suc¬ 


cessful  operation  of  at  least  one  element  of  two  or  noro 


olemonts  in  aotivo  redundancy#  It  is  also  assumed  that  the 


failure  of  any  element  or  elements  in  the  redundancy  places 


no  additional  stress  on  the  remaining  elements.  With  the 


restrictions  staked  above,  the  system  reliability  for  tine 


tQ  can  be  irritten 


W  =  1  -  n  (1-Ri(t0)  ), 


where  1-H^(tQ)  is  the  unreliability  of  the  i™1  oloment  of  m 
in  active  redundancy#  If  the  object  is  to  maximise  system 


reliability  for  time  t  ,  then  it  is  neoossary  to  minimize 


fche  quantity  (1.n(t)) 
i=1  1  0 


for  the  system  elements. 


Since  for  V/eibull  distributed  elements  the  unreliability  can 


be  expressed  as 


1  -  e 


the  object  is  to  minimize 


n  1 

n  (1-  e 

i“1 


V, 


where,  as  in  tho  last  chap tor,  tho  total  expenditure  con- 


iiq#  2k 
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straint  is 

r n  a 


2  A,a,LP.,  1  =  h 

i=1  1  1  1 


If  Equations  21|  and  25  are  used  to  form  tho  Lagrangian 
function,  the  equations  which  result  from  partial  differen¬ 
tiation  are  simultaneous,  non-linear  in  the  GU.  These 
equations  are  not  ameable  to  closed  form  solution.  One 
possibility  in  reducing  tho  complexity  of  these  equations  is 
to  make  the  first  order  approximation 


P 


P 


exp  -to/a  =  i-  to/a  . 

For  element  roliabilitios  greater  than  .72,  this  results  in 
less  than  5  *  error.  With  this  approximation  the  tern  to 
minimize  for  the  system  is 

n  .Pi 
n  fc0 

i=i~37"  • 

Forming  the  Lagrange  function  as 

m  ^  fm  a*  lO^b*  1 

n  *o  -  X  2  a*  a.1/?.  -  r  , 

i=i~a7  Li-i  1  1  1  J 

the  equations  to  satisfy  are 

i  *  A  x  a4-1  n-b, 

-  i  n  _  - Xa,A, a.1  R,  1  =0 

ai  i=i”a7  %  1  1  1 

n  Pi  \  ^  A“bi  “1 

ln(fc0)  n  t0  +  PL  =0 

i=1~a7 


iK i , 2 , • • ,n 


i~1 ,2,#.,n 


“  2  ^a-iPi  -  h  -  o 

i=i  1  1  1 


m(|s  .  26 


Eqs.  27 


Sq.  26 


Combining  Equations  26  and  27  yiolds 


ft- 


iso.  a 


lixnlV 


which  is  a  familiar  result,  Summing  Equation  26  over  all  m 
elements  gives 

331  ai  /V^i 

with  n  a,  a,  & 

i=1  1  x  1 

replaced  by  H,  the  multiplier A  is 

•  a  i] 


This  multiplier  in  Equation  2?  yiolds: 

in(t  )n  -fn  "fill  1 1  b  a6!/?"15!"1 

°  i=1  a±  Li=1  aiJ  biAiai 


Rearranging  in  terms  of  CL^, 


Oj  .  [s  '“(‘olft1 


,VH  1 

ai 


bi  Ai  X  - 

i=i  ai 


Since  Pi  =  bj/a^^  ln(t0)  , 

r  b, 

a,  =  11  <v 

^  b, +1  bi  R  « 

Ai&i  1  ln(t0)  1  Z  i 
1  1  0  1=1  *3 


3q.  sr 

A? 


which  is  the  optimal  CL  allocation  for  each  element  when  m 
different  elements  comprise  the  system. 


:aAHP?&  S 

A  design  is  being  considered  for  a  programmable  elec¬ 
tronic  calculator.  Because  a  problem  can  be  programmed  or 
inputed  via  a  manual  keyboard,  the  programmer  and  keyboard 
are  considered  to  bo  parallel  elements.  Either  clement  can 
fail  and  a  problem  can  still  be  processed  using  the  other 
element  as  an  input  device.  The  cost  function  for  the 
electronic  programmer  has  been  estimated  as 


C 


1 


5.o  .25  0-i.55 
e  a1  ^ 


♦ 


The  cost  function  for  the  keyboard  has  similarly  been  esti¬ 
mated  a3 


If  v2,000  is  to  be  spent  on  the  input  system,  what  0-  and  (2 
should  bo  allocated  to  each  element  to  naxinizo  reliability 
for  5C0  hours  operation? 


Using  Equation  26  with  the  above  oost  parameter*,  the 
optimum  (2 ^  for  elements  one  and  two  are 

/£,  =1.00  ana  /?2  *  2.1(2  . 


iTne  optinun  O..  can  be  found  by  Equation  29  and  are 


a1  =  e 


6.91 


and 


CLg  **  0 


22.65 


These  CL  and  (2 substituted  in  theip  reapootive  oost  functions 
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yield 


^  =  ’*»1  ,070  and  =  ,;930. 


The  system  unreliability  is  m  jj. 

n(tV  a,), 

i=i  0  1 


,™<V> 76  for 


500  hours  operation. 


MI'IICTM  COST  ?0a  PLIABILITY  asqpi.rc-UEra 

L’ho  object  of  the  following  analysis  is  to  minimize 
1  a.  ^-b. 


twp;1 


when  a  system  reliability  requirement  must  be  not.  If  the 

system  reliability  requirement  is  R(t_)  =1-  K(t_)  for  t 

so  o  o 

hours  operation,  the  object  is  to  minimise 


m  a.  ^-b. 

2  a.  a.1/?.  1 

i=i  11  ' 


subject  to 


m 


n  ^o_  =  K(t0). 
i=i  cq  0 


For  this  case  the  La^ran^ian  is 

n  a,  n-b,  rn  .$i  -) 

‘  -  X[n  J.  -KU0jj 


so  that  the  equations  to  satisfy  are 

a,-1  X  m 

Hhai  Pt  +  4-  n  =  0  i=1»2,..,m  , 

1  1  ai  i=i  a* 

a,  rt-b,-1  .  n  ^i 

biAi ®i  t  Xln(t0)II  tG  s  0  i=1*2,..,n  , 

i=1  au 


Ecs.  ; 

30 


Eqs.  ; 

31 


and  m 


n  wo  -  K(tJ  «  0. 


i*i  at 

As  before,  combining  Liquations  30  and  31  results  in 


n  -  bj 

&l  ln(tQ; 


i=1 ,2,.*,n  . 


Substituting  the  reliability  constraint  in  nuation  30  and 


multiplying  by 


gives 


/ai  r_  .  1/ai 


*  ' '  “1  r  i 

’o  _  Ml  1  a.  A,  t  1  1 

ai  '  - - 7K 

L  x<t0)^±lJ 


Taking  the  product  of  this  expression  over  all  m  elononts 


yields 


aiftl  i 


15  1  r  aiKo  i 

■  K(t0)  =  ft  (-  i)ai  *1 

i_1  ai  1=1  I-  K(fc0)>Gii  - 


2  i 

rxii=1i-  i 
rAJ  -  in 


r  b±+1  b,  b.  1 

1  m  aA  Aj^e  ^lnft^)  1  Sj 

w2i - ri - 

°  i1L  KCt  )  bLX  -1 


The  optimum  multiplier  is  thus 


_  rbi+1  b-s  b.-,l"l  9""1 

-X.-iprn  “i1  v  'to'V»  *> 

,vlVi*i  SI - 

L  L  K{tc)  bf  -I  J 

•ith  this  result  tho  o-.tinum  CL  for  the  ith  cl  orient  can  be 


found  as 


2q.  8 


tfq.  33 


hk 


where  refers  to  the  time  for  which  the  system  reliability 
has  been  specified. 


EXAMPLE  6 

Por  the  programmable  calculator  of  Example  what  min¬ 
imum  cost  could  be  expected  to  achieve  a  ,99  reliability 
for  £00  hours  operation?  What  CL  and  (3  should  be  allocated 
the  two  elements? 

The  optimum  level  for  tho  can  again  be  found  from 
Equation  26  and  are 

p,  =  1 .00  and  Pz  -  2.1,2  . 

With  a  K(500)  value  of  .01,  the  optimum  multi  ilior  -X  can 
be  found  in  Equation  3 2  and  is  exp  9o75.  The  resultant 
allocation  for  Ct^  and  CL  ^  are  expl5#42)  and  ex>>(20,38)  from 
Equation  33*  This  results  in  an  unreliability  of  ,00lj.6  for 
element  cr.v  and  2,22  for  element  two.  In  the  development  of 
the  allocation  for  the  parallel  system  no  restriction  was 
imposed  on  the  value  of  unreliability  an  element  could  be 
allocated,  however  a  value  of  unreliability  greater  than  one 
is  physically  meaningless.  When  this  occurn  a  close  look 
should  be  taken  at  the  elements  envolvod  to  see  whether  or 
not  they  are  actually  required  in  the  systom.  If  sof  the 
lowest  reliability  pormissable  should  be  allocated  those  ' 

k 

* 
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elements  in  order  to  maximize  overall  system  reliability* 

For  the  example  under  discussion  it  is  docidod  that  *60  is 
the  lowest  admissible  reliability  for  olement  two;  it  is 
necessary  for  element  one  to  have  .95  reliability  at  500 
hours  to  satisfy  the  system  requirement#  V/ith  the 
hold  constant  the  new  Gj.  are 

7.82  16.00 

G^  »  e  and  G^  =  e 

These  G  and  (2  result  in  element  costs  of  $1,030  for  element 
one  and  $230  for  element  two  for  a  system  cost  of  $1,260# 

For  the  situation  where  more  than  two  elements  comprise 
the  system  and  one  is  allocated  an  unadmis sable  level  of 
unreliability  the  analysis  would  proceed  an  follows*  The 
elements  which  were  allocated  too  large  on  unreliability 
v.’ould  be  set  to  their  minimum  reliability.  The  new  system 
reliability  requirement  could  be  found  by  factoring  out  the 
unreliability  due  to  the  elements  set  at  their  minimum  re¬ 
liability.  A  new  allocation  can  now  be  made  ovor  the  re¬ 
maining  elements. 

OPTIMAL  REDUNDANCY  FOR  IDENTICAL  ELKn^TS 
The  situation  is  often  o.icountered  where  identical  ele¬ 
ments  are  arranged  in  active  redundancy  to  improve  the 
system  reliability.  When  this  is  the  c^so  it  might  be  ad¬ 
vantageous  to  know  how  arny  elomonts  should  be  U30d  t.o 
achieve  the  most  reliable  system  for  a  fixed  expenditure  H. 
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If  a  cost  function  of  the  form  AQa/?“b  is  known  for 
these  elements |  it  is  possible  to  find  the  optimum  redundan¬ 
cy  m*  For  m  elements  comprising  the  systom,  tho  unrelia¬ 


bility 


m  .  A- 
XI  to 


i*i  at 


reduces  to 


!?] 


for  identical  elements*  The  cost  associated  with  this 


system  is 


a/0"b 
mACX  jj 


Thus  for  this  case  the  Lagrangian  may  bo  formed  as 

-  X|nAO. */3-»  -  !l )  . 

The  simultaneous  equations  to  consider  arc 


[$■ 

+  Xamna’'#' 


m  1 

a  L  a  J 


n  ln(tQ) 


«  0  , 


=  0  , 


m 


*  bnua a/3'b_1  =0 


a  ^-b 

nA  CL  f*5  -  H  ®  0  • 


2q.  9 

3* 

2q.  * 
3S- 

2q.  % 

3C 


Simultaneous  solution  of  Equations  3!?  and  36  in  terms  of  f 3 


yiolds 


4q.  37 


liquations  34  and  35  can  be  combined  to  civo 

*tft  -  i 

or  a  13  b+1 

a  *  et 


Sq.  38 


Using  the  cost  constraint  with  the  optimum  CL  and  /? ,  the 


expression  for  the  optimum  redundancy  can  bo  written  as 

b 

II  b 

m  -  s+m - b  • 


2q*  38 


Ae  a  (  ln(tQ)  ) 


Inhere  is  no  guarantee  that  m  will  assume  an  integer  value  in 


Squation  39  although  anything  but  an  integer  number  of  re¬ 


dundant  elements  has  no  physical  significance*  Unloss  the 


optimum  m  is  found  to  bo  and  integer,  the  integer  Just  small¬ 


er  and  just  larger  should  be  tested  to  see  which  is  the  true 


optimum.  This  is  accomplished  by  dividing  the  total  expen¬ 


diture  into  equal  parts  for  each  element  at  each  integer 


tested.  Using  Equation  37  and  the  cost  equation  for  each 
element  the  optimal  01  and  j3  can  be  found.  With  those  two 


parameters  known  for  each  element  the  system  reliability  can 


be  determined  and  the  largest  one  chosen  for  the  integers 


tested. 


2XAHPL3 


A  multi-engine,  all  terrain  vehicle  is  being  considered 
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where  the  estimated  cost  of  achieving  Woibull  parameters  Cl 
and  ft  for  each  power  train  is  expressed  as 

8.5  .Sw-IO.O 

C  =  e  CL  ft  • 

For  preliminary  analysis,  one  power  train  operating  success¬ 
fully  is  considered  adequate  for  mission  success#  If  it  is 
decided  that  <$20,000  should  be  spent  on  the  power  train(s), 
what  is  the  optimum  number  of  power  trains  and  what  power 
train  reliability  can  be  expected  for  1000  hours  operation? 

The  optimum  CL  ,  ft  and  m  can  be  found  from  Equations  37 , 
38  and  39  with  the  cost  parameters  given  above#  The  re¬ 
quired  values  aro 

22.0 

=  2.90  Cl=  0  m  «  2.72  . 

It  is  necessary  to  test  integers  on  both  sides  of  m  to  see 
which  yields  the  lowest  value  of  unreliability.  With  m 
equal  to  2,  each  element  has  a  reliability  of  #933  for  a 
system  reliability  of  .9955*  With  m  equal  to  the  element 
reliability  is  .8l4.fi  for  a  system  reliability  of  .9964# 
Therefore  three  power  trains  is  the  optimal  level  of  redun¬ 
dancy  from  a  cost  standpoint.  Figure  2  illustrates  tho  var¬ 
iation  in  system  reliability  with  the  number  of  redundant 
elements  when  the  total  expenditure -is  $20,000. 
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CHAPTER  V 


SUMMARY  AND  CONCLUSIONS 


A  cost  based  optimum  reliability  allocation  toohnique 
has  been  developed.  Two  basic  configurations  have  been 
considered  for  allocation,  series  systems  and  parallel 
systems*  The  technique  is  based  on  the  minimization  of 
system  unreliability  subject  to  a  total  expenditure  to 
achieve  system  reliability  improvement*  The  dual  of  this 
problem  is  also  considered,  namely  the  minimum  expenditure 
necessary  to  achieve  a  predetermined  system  reliability 
requirement.  The  applicability  of  this  technique  is  con¬ 
tingent  on  the  accuracy  with  which  each  elements  probabil¬ 
ity  density  function  can  be  approximated  by  the  two  para¬ 
meter  V/eibull  distribution.  In  addition  it  is  necessary 

m 

that  parameters  can  be  found  for  the  prupo^  w., ;;  function 
which  discribes  the  cost  associated  with  achieving  a  reli¬ 
ability  level  in  an  olemont*  A  method  is  outlined  in  the 
Appendix  for  estimating  parameters  for  tho  cost  function. 
The  solution  to  the  reliability  allocation  problem  is  found 
by  the  application  of  tho  Lagrange  multiplier  method. 

For  the  case  of  the  series  configuration,  two  snecific 
systems  are  treated.  The  first  is  tho  system  where  each 
element  has  an  exponential  distribution  of  timo  to  failure 
(constant  failure  rate),  and  the  second  13  whore  the  p.d.f. 
of  time  to  failure  for  each  element  oan  be  adequately 
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described  by  the  Weibull  distribution#  Where  the  elements 
have  a  constant  failure  rate,  an  optimum  failure  rate  allo¬ 
cation  is  shown  whore  the  total  expenditure  is  specified  or 
tho  system  failure  rate  is  specified#  For  the  case  where 
all  the  elements  have  a  Weibull  failure  rate  a  closed  form 
solution  is  not  possible  with  the  proposed  cost  function# 

An  iterative  technique  is  required  to  find  a  suitable  multi¬ 
plier  and  thus  the  optimum  allocation#  A  modification  of 
the  cost  function  is  presented  which  allows  tho  explicit 
solution  of  optimum  Weibull  parameters  for  each  element  and 
thus  tho  money  to  be  spent  for  tho  improvement  of  each  ele¬ 
ment.  To  illustrate  the  use  of  the  allocation  technique  a 
numerical  example  is  shorn  for  each  section.  Since  the 
exponential  p.d.f.  is  a  special  case  of  the  Weibull  p.d.f., 
either  approach  can  bo  employed  to  allocate  to  a  constant 
failure  rate,  series  system.  The  choice  should  denend  on 
the  accuracy  of  the  cost  function  in  defining  true  cost 
versus  element  failure  rate. 


For  the  parallel  system  a  method  is  developed  to 
allocrte  Weibull  parameters  to  each  element  to  minimize 
system  unreliability  for  a  given  expenditure.  As  with  the 
series  system  it  is  assumed  that  the  p.d.f.  of  each  element 


can  bo  adequately  described  by  the  Weibull  distribution. 


The  cost  function  is  tho  same  as  used  with  the  previous 
series  model.  The  dual  of  the  above  problem  is  also  con¬ 
sidered  where  a  system  reliability  requirement  must  be  net. 


One  additional  situation  is  considered;  finding  the  optimum 
number  of  identical  elements  in  parallel  to  maximize  system 
reliability  with  an  imposed  expenditure.  To  allow  an  analy¬ 
tical  solution  of  the  allocation  problem  a  first  order  ap¬ 
proximation  is  made  for  element  reliability.  This  approxi¬ 
mation  results  in  less  than  error  for  element  reliabili¬ 
ties  greater  than  .72#  Following  each  development  is  a 
numerical  example. 

One  general  conclusion  that  can  be  made  rogarding  reli¬ 
ability  allocation  is  that  for  anything  but  very  simple 
system  configurations  the  mathematics  involved  becomes  quite 
complicated.  Systems  of  any  complexity  (series-parallel, 
dependent  element)  in  most  cases  require  a  dynamic  program¬ 
ming  approach  and  subsequent  digital  computer  solution. 
Reliability  allocation  is  an  important  step  in  the  concept 
tual,  development  and  design  phases  of  a  system.  As  such  it 
is  a  task  which  must  be  performed.  If  a  model  is  available 
which  can  be  used  to  treat  the  specific  system  under  consid¬ 
eration,  it  should  be  tested. 

SUGGESTIONS  FOR  FURTHER  RESEARCH 

Work  on  the  development  of  the  models  presented  in  this 
paper  could  be  expanded  in  several  directions. 

In  the  series  model  using  the  Weibull  approximation,  it 
iras  found  that  explicit  analytical  solution  was  not  possible 
for  the  optimum  Weibull  parameters.  An  iterative  technique 
is  required  to  locate  the  optimum  Lagrange  multiplier.  A 


digital  computer  program  would  enharoe  the  usefulness  of  the 
allocation  technique  and  serve  as  a  vool  for  sensitivity 
analysis  and  for  study  of  alternative  expenditure  policies. 
Solution  is  not  limited  to  the  digital  computer.  Tno  simul¬ 
taneous  equations  resulting  from  the  Lagrangian  function 
could  be  solved  by  the  use  of  an  analog  computer  and  sensi¬ 
tivity  of  the  optimum  allocation  could  bo  studied  by  the 
variation  of  problem  parameters. 

The  optimum  allocation  for  the  parallel  system  was 
found  only  after  a  first  order  approximation  was  mado  on 
element  reliabilitios.  Perhaps  an  analytical  solution  exists 
for  this  problem  which  does  not  necessitate  this  simplifying 
assumption.  If  not,  it  would  be  worth  while  to  test  the 
closeness  of  the  allocation  developed  in  this  paper  to  the 
true  optimum.  This  could  be  done  by  solving  the  original 
non-linear,  simultaneous  equations  of  the  Lagrangian  by  an 
iterative  technique.  Again,  a  digital  computer  could  be 
employed  to  solve  these  original  simultaneous  equations. 

One  further  area  of  research  is  in  defining  suitable 
parameters  for  the  proposed  co3t  function.  The  method 
presented  in  the  Appendix  is  very  limited  in  that  it  is 
necessary  to  assume  a  fixed  ratio  between  elements  moan  life 
and  standard  deviation  of  time  to  failure.  This  Is  equiva¬ 
lent  to  holding  the  beta  or  shape  parameter  in  the  Ueibull 
distribution  constant.  This  might  not  be  a  realistic  cost 
versus  reliability  assumption  for  certain  elements.  A 


valuable  contribution  to  reliability  work  would  be  the  re¬ 
search  and  cataloging  of  distributional  fonns  of  various 
system  elements  and  components,  and  the  cost  associated  with 
the  reliability  levels  achieved.  Such  a  listing  would 
greatly  aid  in  the  development  of  more  accurate  cost  func¬ 
tions  for  reliability  allocation. 
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APPENDIX 


ESTIMATION  OP  V/EIBULL  PARAMETERS 

The  accuracy  of  the  optimum  allocation  models  that  were 
presented  in  Chapters  III  and  IV  are  strongly  dependent  on 
two  assumptions.  The  first  is  the  assumption  that  the 
p.d.f .*s  of  the  system  elements  can  be  adoquatoly  described 
by  the  V/eibull  p.d.f .  If  a  Weibull  distribution  can  bo  found 
it  is  still  necessary  to  define  a  cost  function  which  re¬ 
lates  the  achievement  of  improved  reliability  to  tho  cost  for 
improvement.  The  key  to  overcoming  both  of  those  problems 
lies  in  the  accurate  estimation  of  the  V/eibull  parameters, 
alpha  and  beta,  for  an  elements  distribution.  Several 
techniques  are  available  for  finding  suitable  V/eibull  para¬ 
meters  to  describe  a  given  distribution.  Three  posr.il  le 
techniques  are  1 )  equating  the  mean  and  variance  of  the 
V/eibull  p.d.f.  to  the  corresponding  parameters  of  the  ele¬ 
ment1®  distribution,  2)  defining  alpha  and  beta  such  that  the 
area  under  the  element’s  p.d.f.  is  equivalent  to  the  area 
under  the  V/eibull  p.d.f.  for  some  time  period  and  3)  use 
maximum  liklihood  estimators  with  points  from  the  elements 
distribution.  Each  of  theso  methods  will  be  discussed  pre¬ 
sently. 

A  very  straightforward  approach  to  obtaining  estimates 
of  alpha  and  beta  i3  to  equato  distribution  means  and  var¬ 
iances.  Por  a  V/eibull  distribution  of  the  form 
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t  /a 


the  mean  tyX  )  and  variance  ( <7  )  can  he  written  as 


-I*2 


and 

1  a  "  1 

These  equations  when  rearranged  in  terms  of  the  tfeibull  para* 
meters  become 


=*  1  + 


(CT/jU.) 


and 


Thus  if  the  ratio  of  standard  deviation  to  mean  is  known, 
beta  oan  be  estimated.  With  beta  and  the  distribution  mean 
the  second  equation  oan  be  used  to  estimate  alpha. 

Another  possibility  for  estimating  an  alpha  and  beta  is 
to  equate  area  under  the  V/eibull  p.d.f.  for  a  given  time  in¬ 
terval  to  the  area  under  the  element's  p.d.f*  If  for  time 
t^  the  system  element  has  reliability  an  equation  of  the 
form 


-t$a 


R, 


oan  be  found  for  all  points  3^ (t^)  along  the  reliability 
curve.  At  least  two  data  points  are  required  to  find  esti¬ 
mates  of  alpha  and  beta.  For  each  pair  of  points  the  two 
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simultaneous  equations 

4°. 


and 


mbp/Q, 

a  *  R2 


age  solved.  The  result  of  this  solution  is 


^  In  (  lnR^lnRg) 
H  In 


and 


a»  -ti 


In?,'" 


If  more  than  two  data  points  are  available  for  parameter 
estimation,  the  points  can  be  taken  in  all  combinations  of 
two  at  a  time  and  equations  solved  for  several  estimates  of 
alpha  and  beta.  For  n  data  points  it  is  possible  to  find 
nl/((n-2) I (2) ! )  estimates  of  the  V/eibull  parameters.  An 
overall  ostimate  of  these  parameters  can  be  made  from  the 

average  of  the  so  obtained  estimates. 

♦ 

Maximum  liklihood  ostlmatore  can  be  used  to  estimate 
alpha  and  beta  if  the  distribution  of  the  element  times  to 
failure  is  initially  assumed  to  be  a  Weibull  p.d.f.  For  a 
sample  of  n  tines  to  failure  (t^)  the  maximum  liklihood 
estimators  for  alpha  and  beta  are  (4) 

nCL-  T  J?  »  0 
i»1  1 

_  n  n 

end  n/£>  +  T  lnt*  -  X  t*  lnt,  *  0  . 

i*1  1  i*1  1  1 
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Solving  these  equations  for  the  V/eibull  parameters  would  of 
oourse  require  an  iterative  procedure.  To  by-pasr,  this 
difficulty  beta  can  be  found  by  solving  the  equations 


F 


and 


n 

n  Cl  *•  2  t.  ~  0 
i=1  1 


n 

2  t 
i=i  1 
n 


simultaneously.  This  second  equation  is  an  unbiasod  esti¬ 
mate  of  the  distribution  mean.  Solution  of  these  equations 
leads  to 

2.  fi 

2  t±  =o. 


l  r=i 


From  this  equation  an  estimate  of  beta  can  be  found,  An 
estimate  of  alpha  is  then  found  from 

a  n 

«•&,* *  • 


The  use  of  one  particular  parameter  estimation  tech¬ 
nique  over  another  is  dependent  on  several  faotors#  If  the 
?,d,f*  is  known  for  the  element  under  question,  the  first 
method  should  be  used,  Thr,  in  an  ard  Variance  of  the  distri 
bution  can  be  equated  and  the  desired  CL  and  (2  parameters 
determined,  !fnen  only  two  points  are  available  on  the  ele¬ 
ment  reliability  curve,  the  second  estimation  method  is 
easily  used  and  accurate  enough  for  preliminary  allocation* 
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tlhen  a  sufficient  volume  of  data  is  available  and  the  V/oi- 


bull  p.d.f.  is  assumed  then  maximum  liklihood  estimators 


are  the  most  aocurate  estimators  available. 


3&TIKATI0N  0?  COST  FUNCTION  PAdAKJSVidS 


jtflM  The  following  is  a  technique  whereby  the  parameters  A, 

\  ■  * 


a  and  b  can  be  estimated  for  the  cost  function 


a/rb 
a  a  p 


One  simplifying  assumption  of  this  technique  is  that  the 


ratio  of  element  moan  life  to  standard  deviation  of  times  to 


failure  is  constant.  This  is  equivalent  to  holding  the 


"shape"  of  the  distribution  constant,  For  a  preliminary 


l  location  this  is  not  considered  a  serious  drawback.  It 


seems  reasonable  that  even  though  the  moan  life  of  an  equip¬ 


ment  increases  the  general  shape  of  its  failure  probability 


distribution  should  not  be  altered  greatly.  The  accuracy  of 


this  simplifying  assumption,  however,  will  bo  born  out  by 


the  Noibuli  parameters  estimated  for  the  element  in  question. 


In  addition,  at  least  two  oost  estimates  are  reouired  for 


each  element  on  the  cost  vorsus  element  reliability  curve. 


If  the  proceeding  requirements  are  met  the  procedure  is  as 


follows.  Using  an  estimating  technique  mentioned  earlier, 
h’eibull  parameters  Q»1#  CL^  and  p  are  found  for  the  element. 
Only  one  /?  is  defined  since  the  ratio  of  mean  life  to  stan¬ 


dard  deviation  remains  oonstant.  Associated  with  theie  two 


reliability  levols,  S(CLj,$,t)  and  R(GL>,/?,t),  for  oper- 


j 

1 


6t 


ating  tine  t  are  two  coat  estimates  C ^  and  C^.  The  cost 
equations  resulting  from  these  points  on  the  reliability 
curve  aro 


a  0-b 

and  C2  ts  xaz(S 

where  A  is  the  position  parameter  and  a  and  b  aro  the  slope 
parameters.  Taking  the  natural  logarithm  of  both  equations 
and  combining  yields 


.  ln(0v/C2) 


Since  it  was  found  for  the  analyses  in  Chapters  III  and  IV 
that  the  optimal  level  of  beta  for  all  elements  17111  bo 
b/a  ln(t0),  b  can  be  estimated  as 
b  =  /?a  ln{t0) , 

whore  as  beforo  tQ  io  the  system  operating  tino  ovor  which 
reliability  is  being  optimized.  Using  an  estimation  tech¬ 
nique  previously  discussed,  CL  and/?  can  bo  doteminod  for 
two  reliability  and  cost  levels  of  the  element,  With  tuo 
Cl  values  defined,  the  slope  parameter  a  can  bo  found  from 
Souatior.  1A.  Once  a  i3  dote  x.oined,  b  can  be  found  from 
Equation  2A.  With  both  a  and  b  determined,  tho  oquution  for 
the  position  paraneter  A  can  be  written  as 


Eq.  * 
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